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SURFACE CURVES FOR STEADY 
NONUNIFORM FLOW 


ROBERT JUN. ASCE 


Few problems the field hydraulic research have been more challenging 
baffling than that concerned with plotting the surface curve assumed 
liquid varied nonuniform flow. Several investigators have sought the 
solution attempting integrate the varied flow equation. This paper, 
however, makes use the differential equation direct application the 
surface slope that defines. The paper confined the study steady flow 
prismatic trapezoidal channels, with the rectangular section considered 
limiting case. Surface profiles plotted this method are compared with 
those obtained integration. 


INTRODUCTION 


quick and simple method determining the water surface profile 
open channels carrying nonuniform flow has been the object search 
engineers and mathematicians for number years. Although many ex- 
haustive efforts have been made solve this problem integration, 
entirely satisfactory answer has been found. The method advocated Boris 
Hon. ASCE, has been applied with more confidence than 
has been accorded the works other researchers, but even this approach has 
its deficiencies that are generally recognized. 

The procedure outlined herein is, regrettably, not intended final 
and complete solution. certainly not presented substitute for the 
scholarly achievements Mr. Bakhmeteff, nor detract from the works 
other earnest investigators. The paper essentially an-enlargement upon the 
obvious relationship defined the differential equation varied flow, and 

comments are invited for publication; the last discussion should submitted 
April 1952. 

Eng. Associate, Bureau Engineering, City Los Angeles, Calif. 
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offered illustration the results obtainable direct plotting the water 
surface slope. 
VARIED 


The flow liquid open channel said uniform when factors 
such bottom slope, depth, velocity, and cross-sectional area remain constant 
from point point. Uniform motion can considered limit that 
approached, for example, long reach channel whose slope just equal 
the rate loss energy. general, therefore, the case that must 
considered the one which the contributing factors vary such manner 
that the energy gradient, the 
bottom slope, 
quently the water surface are 
not parallel. 

The following treatment 
devoted steady nonuni- 
form (or varied) flow pris- 
matic channels. The slope 
either constant varying 
gradually that the effect 
centrifugal action can neg- 
lected. That is, the pressure 
distribution taken hydro- 
static. the calculation 
energy losses, application 
made the Manning version 

entails the assumption that 

the energy gradient varied flow parallel that corresponding uniform 
flow under the same conditions depth, roughness, and velocity. 

the distance along the channel parallel the bottom for ele- 
mentary length channel (dz), illustrated Fig. 1(a), the total head (h) 
any point the stream 


Horizontal Datum 


(a) 


assuming the velocity (V) constant with respect depth. Eq. 
the elevation the channel bottom above horizontal datum, the 
depth water the channel, and the acceleration gravity. 


which coefficient; the slope the energy gradient; and the 


Fak 
oy ‘ 
(1) 
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hydraulic radius. Also, from Fig. apparent that which 


the slope the channel bottom. Referring Fig. infinitesimals 
higher than the first order are neglected, dy, which the surface 
width and the cross-sectional area the channel. 


(3) 


Substituting those derived terms into Eq. gives 


transposing terms, the following final expression for the slope the 
water surface obtained: 


Now, if, indicated Fig. the tangent the channel side slope 
the bottom width the channel, and the surface width 
Also, the wetted perimeter 
The expression for the hydraulic radius therefore 


(6) 


For prismatic channel, constant. the Manning formula, 
1.486 2.208 


inserting the foregoing factors, Eq. can now written 


2.208 
7] n 
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or, 


The constant terms Eq. can denoted 


the indicated substitutions, the differential equation can changed 


the form 


that, for this shape, 


Since for prismatic channel the terms and not vary, the expres- 
sions within brackets Eq. are functions the depth alone. These 
variables the numerator and denominator are designated and respec- 
tively. 


(9) 


and 


The relation for the trapezoidal section can now restated 


The depth functions, and have been evaluated and tabulated for common 
cases side 
PROCEDURE 
illustrate the direct use of.the surface slope solving problems 
varied flow, the following elementary examples are given: 
for the values the factors and are available for the cases vertical side walls and 


side walls having slope and 1.5 with bottom widths ranging from 5-ft inter- 
vals and depths from 0.5-ft intervals. The cost these tables $1.00. 
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Example 1.—Water flowing from under sluice into rectangular canal 
whose bottom width ft. The channel grade 0.0010 and the value 
Manning’s 0.013. The sluice regulated discharge 1,080 per sec 
with depth the vena 


contracta 1.5 ft. The TABLE SURFACE 
The critical depth 
found 4.49 ft. Flow Depth Factor 
under the conditions stated 
will usually lead hy- 2.0 +0.00264 
draulic jump, starting be- 
tween 1.5 and critical 57.30 58.31 +0.00371 


depth 4.49. For the 
purpose illustration, how- 
the curve will computed the assumption that the tailwater low 
enough prevent occurrence the jump. 


2 2 2 


2.208 
and 11, the values and can found. Solution Eq. then 
completed. 
5.0 


c 

3.0 

Bakhmeteff 

a 

a 


300 400 500 
Distance Along Channel Feet 


should remembered that the water slope with respect the 


channel bottom. conventions the calculus, its sign positive with 
increasing depth the direction flow. 


4.0 
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the manner described above, values for intermediate depths have 


been determined and are listed Table 
The water surface curve described these slopes plotted Fig. 
order attain uniform results the sketching curves this paper the 
tangents have been made 


TABLE SURFACE intersect the mid-points 
the depth increments, that 
1481 nat part the curve where 
Varying rapidly. 


For comparison, the sur- 
face plot obtained the Bakhmeteff method shown the dashed curve 
Fig. 

Example experiment conducted the University Illinois 
Engineering Experiment Urbana, 40.39 per sec water 
was caused flow through rectangular channel wide. The flume was 
laid slope 0.0030 and was artificially roughened obtain value 
backwater gate installed the lower end the reach was 


Depth Water (y) Feet 


0.0000 


1600 1400 1000 


regulated this particular test that the depth the downstream gage 
was 4.0 ft. The normal depth for the given conditions found 
equal 1.84 ft. Also, using procedures already described, 0.267 and 
50.7. 


Lansford and Mitchell, Bulletin Series No. 381, Univ. Illinois, Urbana, 1949. 


= 
4 g\ve? 
‘ 


NONUNIFORM FLOW 


The slope calculations are summarized Table 
Fig. the surface plotted means these tangents shown com- 
pared with the curve observed the laboratory. 


CoNCLUSIONS 


Direct application the water surface slope defined the varied flow 
equation may expected yield results close agreement with those ob- 
tained accepted integration procedures and trial-and-error methods. 
Any these approaches involves uncertain amount inaccuracy due 
the assumptions hydrostatic pressure and nondivergent motion. The 
error introduced through determination energy loss the uniform-flow 
formula probably not appreciable except relatively steep sharply non- 
parallel movement. 

Profiles plotted any procedure are inaccurate the vicinity the con- 
trol section for most cases, since the vicinity the critical depth the effect 
centrifugal forces usually considerable importance. The water slope 
this point, when channel conditions permit, can approximated curve 
sketched between adjacent sections upstream and downstream. 
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